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ABSTRACT 

Gas-phase complex organic molecules have been detected toward a range of high- and low-mass star- 
forming regions at abundances which cannot be explained by any known gas-phase chemistry. Recent 
laboratory experiments show that UV irradiation of CHsOH-rich ices may be an important mechanism 
for producing complex molecules and releasing them into the gas-phase. To test this ice formation 
scenario we mapped the Bl-b dust core and nearby protostar in CH3OH gas using the IRAM 30m 
telescope to identify locations of efficient non-thermal ice desorption. We find three CH3OH abundance 



peaks tracing two outflows and a quiescent region on the side of the core facing the protostar. The 
CH3OH gas has a rotational temperature of ~10 K at all locations. The quiescent CH3OH abundance 
peak and one outflow position were searched for complex molecules. Narrow, 0.6-0.8 km s _1 wide, 
HCOOCH3 and CH3CHO lines originating in cold gas are clearly detected, CH3OCH3 is tentatively 
detected and C2H5OH and HOCH2CHO arc undetected toward the quiescent core, while no complex 
molecular lines were found toward the outflow. The core abundances with respect to CH3OH arc 
-2.3% and 1.1% for HCOOCH3 and CH 3 CHO, respectively, and the upper limits are 0.7-1.1%, which 
is similar to most other low-mass sources. The observed complex molecule characteristics toward 
Bl-b and the pre-dominance of HCO-bearing species suggest a cold ice (below 25 K, the sublimation 
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temperature of CO) formation pathway followed by non-thermal desorption through e.g. UV photons 
traveling through outflow cavities. The observed complex gas composition together with the lack 
of any evidence of warm gas-phase chemistry provide clear evidence of efficient complex molecule 
formation in cold interstellar ices. 
Subject headings: astrochemistry — ISM: abundances — ISM: molecules — astrobiology 

1. INTRODUCTION 

Complex organic molecules have been detected toward a range of astrophysical environments, including low-mass 
protostars (van Dishoeck et al. 1995; Cazaux et al. 2003; Bottinclli et al. 2007); however, the origins of these complex 
molecules as well as their fates are uncertain. Commonly-suggested formation routes for the detected molecules include 
various gas-phase reactions starting with thermally evaporated CH 3 OH ice, atom-addition reactions on dust grains, 
and UV- and cosmic ray-induced chemistry in the icy grain mantles that form during the pre-stellar stages (Charnlcy 
et al. 1992; Nomura & Millar 2004; Herbst & van Dishoeck 2009, for a review). The focus is currently on an ice 
formation pathway (e.g. Garrod & Herbst 2006; Garrod et al. 2008) because of the failures of gas phase chemistry 
to explain the observed abundances of some of the most common complex molecules toward low-mass protostars, 
especially HCOOCH 3 . 

Recent experiments on the photochemistry of CHaOH-rich ices have shown that 1) UV irradiation of CH3OH ices 
at 20-70 K result in the production of large amounts of the complex molecules observed around protostars and 2) the 
chemistry has a product branching ratio which is temperature and ice composition dependent (Oberg et al. 2009b). A 
key result is that HCOOCH3 and other HCO-X ices are only abundantly produced in CO:CH30H mixtures. A similar 
result has been reported for proton-bombarded CH3OH and CtLjOILCO ices (Bennett & Kaiser 2007). CO:CH30H 
ice mixtures are probably common since CH3OH forms from hydrogenation of CO (e.g. Watanabe et al. 2003; Cuppen 
et al. 2009). CO evaporates at 17-25 K on astrophysical time scales (Bisschop et al. 2006), and thus HCO-X ices will 
mainly form in cold, UV-exposed ices. At higher temperatures, closer to the protostar, UV irradiation of the remaining 
pure CH3OH ice will instead favor the production of C2H5OH, CH3OCH3 and (CH20H)2- 

In light of this proposed formation scenario, protostars rich in CH3OH ice are natural targets when searching for 
complex molecule sources. The low-mass protostar Bl-b is such a source. From previous SCUBA and 3 mm continuum 
maps, Bl-b consists of two cores Bl-bN and Bl-bS, separated by 20" (Hirano et al. 1999). The cores are similar with 
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Tdust ~ 18 K, M — 1.6 — 1.8 Mq and Lboi = 2.6 — 3.1 Lq. The Spitzer Space Telescope only observed one protostar, 
to the south-west of the Bl-bS dust core at 03:33:20.34, +31:07:21.4 (J2000), but it was still named Bl-b - to avoid 
confusion it will be referred to as the 'protostar' position. From the Spitzer c2d (From Molecular Cores to Planet 
Forming Disks) ice survey the Bl-b protostar has a CH3OH ice abundance of 11% with respect to H2O ice (Boogert 
ct al. 2008), corresponding to ~ 5 x 10~ 6 with respect to H2. The fractional abundance of H 13 CO + decreases by a 
factor of 5 toward the SCUBA core, indicative of CO freeze-out (Hirano ct al. 1999). A HCO-dominated complex 
chemistry is thus expected where the ice is exposed to UV radiation. 

The Bl-b region is complicated by a number of outflows, which may enhance the UV field through shocks and open 
cavities through which stellar UV photons can escape (J0rgensen et al. 2006; Walawender et al. 2009; Hiramatsu et al. 
2010). The outflow of most interest to this study runs in the south-west direction from Bl-b protostar, though it is 
not known whether Bl-b is its source or just happens to lie in its path. The outflow is not observed north-east of the 
Bl-b protostar, where the Bl-b SCUBA core is situated. It is therefore unclear whether the outflow terminates at the 
protostar or continues into the core, hidden from view. If it does penetrate into the SCUBA core, this may increase 
the UV flux in the region by orders of magnitude, enhancing both the UV photochemistry and photodesorption of 
ices. 

The latter is important, since once formed, the complex molecules must be (partly) released into the gas phase to 
be observable at millimeter wavelengths. The second reason for targeting Bl-b is the detection of a large CH3OH gas 
phase column density of ~2.3-2.5x 10 14 cm" 2 toward the protostar (Oberg et al. 2009a; Hiramatsu et al. 2010) and thus 
ice evaporation. The narrow line width and low excitation temperature of the observed CH3OH suggest non-thermal 
ice evaporation - the abundances are consistent with UV photodesorption. Assuming that complex ice chemistry 
products desorb non-thermally at an efficiency similar to CH3OH (^10~ 3 molecules per incident UV photon (Oberg 
et al. 2009b)), the protostellar environment should contain a cold gas-phase fingerprint of the complex ice composition; 
Bl-b may offer the first possibility to observe the earliest stages of complex molecule formation, untainted by either 
warm ice chemistry or gas-phase processing. 

Addressing these predictions, we present IRAM 30-m telescope observations of a 80" x 80" CH3OH map of the Bl-b 
protostellar envelope and nearby dust cores, followed by a search for HCOOCH3, CH3CHO, CH3OCH3, C2H5OH and 
HOCH2CHO toward two of the identified CH3OH peaks. The CH3OH gas abundances and the complex abundances 
and upper limits with respect to CH3OH are then discussed in terms of different ice desorption mechanisms and 
complex molecule formation scenarios. The results are finally compared with complex molecule observations toward 
other low-mass protostars and outflows. 

2. OBSERVATIONS AND DATA REDUCTION 
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The Bl-b protostellar envelope and the nearby SCUBA cores were mapped in three CH3OH transitions in February 
2009 using the HEterodyne Receiver Array HERA (Schuster et al. 2004) on the Institut de RadioAstronomie Mil- 
limetrique (IRAM) 30-m telescope together with the VESPA autocorrclator backend (Tables 1 and 2). The 80" x 80" 
map is centered on R.A 03:33:21.3 and Dec 31:07:36.7 (J2000), the estimated midpoint between Blb-S and Blb-N 
(Hirano et al. 1999). VESPA was used with a 320 kHz channel spacing resulting in a velocity resolution of 0.78 km s _1 
(channel spacing <5V of ^0.39 km s _1 ) at 242 GHz. Observations were carried out through two separate undersamplcd 
(8" spacing) raster maps using beam-switching with a throw of 200" in azimuth. The two maps were off-set by +4" , + 
4" and the combination of the two maps thus results in a single over-Nyquist-sampled map - the half power bcamwidth 
of the 30-m telescope at 242 GHz is ~10" . The rms is ~50 mK in the first map and ^130 mK in the second map in a 
0.39 km s _1 channel. The average rms over the entire map was estimated to ~90 mK. The system temperature varied 
between 600 and 1200 K. 

The second set of observations aimed at detecting complex molecules were carried out in June 2009 with the IRAM 
30-m telescope, using the new EMIR receiver (Table 1). The positions used for pointing correspond to the CH3OH 
maximum toward the center of the SCUBA core at 03:33:20.80, 31:07:40.0 (J2000), termed the 'core' position, and 
a local maximum at 03:33:21.90, 31:07:22.0 (J2000) termed the 'outflow' position. The targeted lines are listed in 
Tabic 3 and the frequencies are taken from the JPL molecular database and the Cologne Database for Molecular 
Spectroscopy (Miiller et al. 2001). The focus is on lines with low excitation energies (15-100 K) because of the 
previously discovered low excitation temperature of CH 3 OH toward the same source. The chosen settings contain at 
least two lines for HCOOCH 3 , CH 3 CHO, CH 3 OCH 3 and C 2 H 5 OH (and additional HOCH 2 CHO and C 2 H 5 CN lines) 
to allow us to constrain the excitation temperatures of any detected species. The observations were carried out using 
three different receiver settings, each including a combination of the EMIR 90 and 150 GHz receivers (Table 1). At 
these wavelengths, the beam sizes are ~27 and 19" , respectively. All three settings were used toward the core position, 
and two toward the outflow position. Each receiver was connected to a unit of the autocorrclator, with a spectral 
resolution of 40 MHz and a bandwidth of 120 MHz, equivalent to an unsmoothed velocity resolution of ~0.1 km s _1 . 
Typical system temperatures were 100-150 K. All observations were carried out using wobbler switching with a 140" 
throw (100" in azimuth and 100" in elevation). 

For both sets of observations, pointing was checked every ^2 hr on J0316+413 with a typical pointing accuracy of 
~2-3" . All intensities reported in this paper are expressed in units of main-beam brightness temperature, which were 
converted from antenna temperatures in paKo using reported main beam and forward efficiencies (B e gf and F e ff) of 
~75 and 95% with the E90 receiver, ~69 and 93% with E150 receiver and a main beam efficiency of 52% at 2 mm for 
the HERA receiver. The rms in mK are reported in Table 1. 

The data were reduced with the CLASS program, part of the GILDAS software package (see 
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http://www.iram.fr/IRAMFR/GILDAS). Linear (first-order) baselines were determined from velocity ranges with- 
out emission features, and then subtracted from the spectra. Some velocity channels showed spikes, which were 
replaced by interpolating the closest two good channels in the map, while no correction was made to the complex 
molecule spectra - there were no spikes close to any expected line positions and thus no correction was necessary. The 
complex organics observation were smoothed to velocity resolutions of 0.4-0.8 km s _1 toward the core position and 2-4 
km km s _1 toward the outflow position respectively to achieve a maximum signal-to- noise while still barely resolving 
the expected lines. 

3. RESULTS 

3.1. CH 3 OH 

Figure 1 shows the integrated CH3OH map superimposed on a SCUBA dust map of the Bl-b region. The dust map 
contains one elongated core, rather than the two sources observed by Hirano et al. (1999). The SCUBA core in the 
remainder of the paper refers to the dust core in Fig. 1, while the 'core' position refers to the nearby CH3OH maximum 
probed for complex molecules. The dust and molecular maps share some features, such as the N-S elongation and the 
'tail' toward the south-west, but the centers of the two maps are offset - the CH 3 OH column density peaks in between 
the SCUBA core and the 'protostar' suggestive of ice desorption on the west side of the dust core due to irradiation 
from the protostar traveling through outflow cavities, as discussed in more detail below. A velocity channel map of the 
same CH3OH line shows that in addition to the quiescent, elongated CH3OH emission that overlaps with the western 
half of the dust emission, there are at least two high velocity components, suggestive of two perpendicular outflows 
which are either directed along the line of sight or which are too young to have any resolved component in the plane 
of the sky (Fig. 2). The two outflow positions are almost coincidental with Bl-bS as reported by Hirano et al. (1999) 
and may thus be evidence for a very young protostellar object embedded in the Bl-bS core. 

Figure 3 shows the extracted spectra at two of the CH3OH intensity peaks in the channel map - the 'core' and 
southeast 'outflow' positions. The CH3OH emission toward the 'core' consists of narrow, symmetric lines with small 
wings. Overall these spectra seem to be dominated by quiescent gas. In contrast the CH3OH lines toward the 'outflow' 
position are more than four times broader and asymmetric, consistent with the assignment to an outflow. The spectra 
toward the 'core' are similar to those observed by Hiramatsu et al. (2010) toward the same line of sight. They did 
however not find any evidence for outflow activity toward Bl-b. The lack of channel maps in Hiramatsu et al. (2010) 
prevents direct comparison, but the detected outflows here are on a much smaller scale compared to the outflows they 
do observe and may thus have been overlooked. 

Table 2 lists the observed CH3OH transitions as well as the line FWHM and integrated intensities toward the 
core and outflow positions. To facilitate comparison with observations of complex molecules and previous CH3OH 
observations at lower frequencies, the CH3OH intensities were extracted both directly from the HERA maps and from 
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maps convolved with 19" and 27" beams, the IRAM 30-m beam sizes at 131 and 89 GHz, respectively. Because of the 
asymmetry in the lines, the integrated intensities were calculated by integrating from -10 to +5 km s _1 around each 
line center, rather than integrating the fitted Gaussians used to estimate the FWHM. 

The excitation temperature of CH3OH toward each position is estimated using rotational diagrams (Goldsmith & 
Langer 1999) in Fig. 4 based on the CH3OH map intensities convolved with a 27" beam. The 'core' and 'outflow' 
positions only contain three CH 3 OH lines each and the derived temperatures are therefore uncertain by at least a 
factor of two. The best fits results in 10 K for the 'core' and 8 K for the 'outflow'. For the 'protostar' position the 
new data are combined with previous IRAM 30-m data observed at ~97 GHz (Oberg et al. 2009a). The resulting 
excitation temperature is 11 ± 1 K, confirming the estimated rotational temperatures toward the 'core' and 'outflow'. 
At the densities expected in molecular cloud cores, CH3OH is easily sub-thermally excited - rotational temperatures 
below 25 K are expected at densities up to 10 6 cm~ 3 even if the kinetic temperature is above 100 K (Bachiller et al. 
1998). NH3 maps of the same region however confirm that in this case the kinetic temperature is close to the CH3OH 
rotational temperature, ~10-12 K (Bachiller et al. 1990). 

3.2. Search for complex molecules 

Building on the CH3OH map, we searched for complex molecules toward the 'core' and 'outflow' positions. All 
acquired spectra are shown in the Appendices. Figure 5 shows six blow-ups focusing on the targeted complex molecules. 
The two lowest lying transitions of HCOOCH3 (E up = 18 K) and of CH 3 CHO (E up = 28 K) are detected toward 
the core region (Fig. 5a, d) at the 5-cr level. Higher lying transitions of the two molecules are not detected at the 3-ct 
level, nor are any transitions of C2H5OH, HOCH2CHO or C2H5CN. One 2-<r transition of CH3OCH3 is tentatively 
seen and its intensity is consistent with non-detections of the other lines. No complex molecule lines are detected 
toward the outflow (Fig. 5). The two securely detected species, HCOOCH3 and CH3CHO, both display narrow line 
profiles between 0.5-0.9 km s _1 . The HCOOCH3 lines show no evidence of asymmetry, while the CH3CHO lines have 
high- velocity wings suggestive of an emission contribution from the '02' outflow. 

All data for the targeted lines are listed in Table 3, including molecule transitions, catalogue frequencies, E\ ow , 
FWHM of detected lines, integrated intensities for detected lines and integrated intensity upper limits for non- 
detections. The upper limits on higher lying HCOOCH3 and CH3CHO lines are consistent with the ~10 K excitation 
temperature of CH3OH toward the same line of sight. Thus a 10 K excitation temperature is assumed when calculating 
column densities. 

Using these excitation temperatures, the total HCOOCH3 column density (A + E) toward the Bl-b core is 8.3 ± 
2.8 x 10 12 cm -2 , averaged over the telescope beam, and the total (A + E) CH3CHO column density is 5.4 ± 1.4 x 10 12 
cm~ 2 , which corresponds to abundances with respect to CH3OH of 2.3 and 1.1% respectively. The different beam sizes 
at different frequencies are accounted for by comparing the HCOOCH3 to the CH3OH column density derived with a 
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27" beam and the CH3CHO to the CH3OH column density derived with the 19" beam. The upper limit abundances 
for the other molecules are <0.8, <1.0% and <1.1% for CH3OCH3 , C 2 H 5 OH and HOCH 2 CHO toward the core. The 
most significant upper limit toward the outflow is the total HCOOCH3 abundance of <1.4% with respect to CH 3 OH. 
The calculated abundances and upper limits of complex molecules arc summarized in Table 5. The lack of detections 
toward the outflow does not prove that it contains significantly less complex molecules than the 'core', but the outflow 
position is definitely not richer in complex organics than the core region. The outflow upper limits are comparable to 
the detected abundances toward another low-mass outflow L1157 (Arce et al. 2008) 

4. DISCUSSION 
4.1. The origin of the CH 3 OH gas 

The CH3OH average abundance toward the quiescent CH3OH 'core' is ~ 2 x 10 -9 with respect to Nh 2 , assuming a 
SCUBA dust intensity conversion factor of N H2 /S£ eam of 1.3xl0 20 cur 2 (mJy beam- 1 )- 1 (Kauffmann et al. 2008). 
This abundance is too high to explain with gas-phase formation of CH 3 OH through any mechanism investigated so 
far, i.e. gas phase reactions in quiescent cores result in CH3OH abundances below 10 -13 (Garrod et al. 2007). The 
observed CH 3 OH gas is thus a product of ice desorption. 

The protostar Bl-b belongs to a small set of young low-mass stellar objects with CH3OH ice abundances of >10% 
with respect to H2O (Boogert et al. 2008, Bottinclli et al. submitted to A&A). CH3OH ice has been previously 
observed to vary on small scales, however, and the ice abundance towards the 'core' and 'outflow' may be as low as 
1% with respect to H2O ice. Still, <0.5% of the CH3OH ice must evaporate to account for the observed CH3OH 
abundance of ^2xl0 -9 toward the 'core', assuming a H2O abundance of ~ 5 x 10~ 5 with respect to H2. 

The extent of the CH3OH gas (~10000x 15000 AU) together with the low excitation temperature of CH3OH and 
the offset between the CH3OH 'core' and the 'protostar' exclude thermal ice evaporation as a major contributor to 
the observed abundances - a low-mass protostar cannot heat more than a few 100 AU to the sublimation temperature 
of CH3OH of ~80 K. This leaves non-thermal evaporation of ices, which can be broadly divided into two categories: 
ice evaporation due to grain sputtering in shocks and less violent non-thermal ice evaporation due to UV irradiation, 
cosmic rays or release of chemical energy (Jones et al. 1996; Shen et al. 2004; Garrod et al. 2007). Grain sputtering is 
often invoked to explain excess CH3OH in protostellar outflows (e.g. Bachiller et al. 2001) and is probably responsible 
for the factor of a few higher abundances of CH3OH gas in the two outflow positions around the dust core, using the 
same conversion factor as above to derive the dust column densities from the SCUBA map. The narrow CH3OH lines 
in the remainder of the map suggest a different origin for the gas toward the protostar and 'core' regions. 

The asymmetric distribution of CH3OH compared to the dust core further suggest that ice is desorbed by irradiation 
from the direction of the protostar, illuminating the south-west side of the core. The exact nature of the irradiation is 
difficult to constrain without detailed modeling, but UV photodesorption is known to be efficient (Oberg et al. 2009b). 
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The dust core has a visual extinction of >100 mag. For radiation to affect the 'core' region thus requires extensive 
cavities in the interstellar medium. As discussed in the Introduction, there are several large-scale outflows crossing the 
Bl region, which could have carved out such cavities (J0rgensen et al. 2006; Walawender et al. 2009). In particular 
Fig. 1 shows that there is a visible outflow in the 'right' direction originating either from the Bf-b protostar or from 
the nearby Bl-d protostar to the southwest. In cither case, the outflow may continue on towards the 'core', hidden 
from observations by the large column of dust, resulting in UV irradiation of the icy grains on the side of the dust 
core closest to the protostar. 

UV irradiation can be generated in a number or ways within such an outflow cavity. If an outflow is still active in 
the direction of the 'core' position, the UV radiation may originate in the jet shock itself (Reipurth & Bally 2001). 
Some outflows are also known to become hot enough to emit x-rays (Pravdo et al. 2001), which may directly desorb 
the ice or desorb it through secondary UV photons. The narrow line width of the CH3OH emission suggest however 
that the CH3OH gas docs not originate close to shocked gas. It seems more likely that a past outflow has opened a 
cavity between the protostar and the 'core' position and that UV radiation from this neighboring low-mass protostar 
protostar is desorbing the ices. Low-mass protostars are known to have excess UV fluxes compared to the interstellar 
radiation field, originating in the boundary layer between the protostar and the accretion disk (Spaans et al. 1995). 

UV photodesorption is thus a probable source of the CH3OH gas distribution towards the Bl-b dust core. Observa- 
tions of UV tracers are however needed to confirm this scenario, since we cannot completely rule out that the CH3OH 
are the left-overs from grain-sputtering by the same shock that opened up the cavity, if it passed through the area 
10 4 -10 5 years ago, the typical depletion time-scale at molecular cloud densities. 

4.2. The origin of complex molecules toward Bl-b 

The similar line widths and excitation temperature upper limits of HCOOCH3 and CH3CHO compared to CH3OH 
suggest that the complex molecules toward the Bl-b core originate from photodesorbed ice as well. This is a new 
potential source of complex organic molecules around low-mass protostars compared to what has been previously 
suggested in the literature, where complex molecules have been observed in a small warm core or disk where ices have 
evaporated thermally (e.g. Bottinclli et al. 2007; J0rgensen et al. 2005) or in shocks following grain sputtering (e.g. 
Arce et al. 2008). In difference to sputtering and thermal desorption, which quickly destroys the entire ice mantle in 
a small region, UV photodesorption can release a small fraction of the ice over large quiescent regions. This offers 
the opportunity to study ice chemistry in situ as it evolves around protostars, as previously suggested in Oberg et al. 
(2009a,b). 

Non-thermal desorption of HCOOCH3 and CH3CHO is only a probable source of the observed gas if the molecules 
can be produced efficiently on grains, however. As discussed in the Introduction, such a production channel exists. UV 
processing of COiCFFjOH-rich ices produces complex molecules at abundances that are consistent with the observed 
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abundances of a few % with respect to CH3OH. 

Moreover this formation path explains the relative abundances of HCOOCH3 and CH3OCH3 or C2H5OH. Figure 6 
illustrates how the addition of CO affects the complex CH 3 OH chemistry to produce the observed molecules. Quanti- 
tatively, UV irradiation of a CH30H:CO 1:10 ice mixture at 20 K under laboratory conditions results in the production 
of a CH3OCH3/HCOOCH3 ratio of <0.1 and irradiation of a pure CH 3 OH ice in CH3OCH3/HCOOCH3 ratios of >1.3. 
While the laboratory results cannot be directly extrapolated to Bl-b, a CO-rich ice is certainly needed to produce the 
observed low CH30CH3(C2H50H)/HCOOCH3 ratios; without CO in the ice, the most recent model predict ratios of 
>10, rather than the observed <0.5 (Garrod et al. 2008). 

4.3. Comparison with previous complex chemistry observations 

With energetic CO:CH30H ice processing followed by non-thermal desorption established as the most probable 
explanation for the observed complex molecules in Bl-b, it is interesting to see how the chemistry toward Bl-b 
compares with complex chemistry observations toward other low-mass sources: NGC1333-IRAS 2A, 4A and 4B, IRAS 
16293-2422 and L1157. All sources have been observed with single-dish telescopes where the emission region is assumed 
to be unresolved. These observations thus contain emission from both hot and cold material - dependent on the density 
profile one or the other component may dominate. In addition IRAS 16293-2422 has been studied extensively using 
intcrfcromctry. In these observations two cores, A and B, are resolved and emission from the cold envelope can be 
excluded. 

Because of assumptions on the emitting region for complex molecules in most previous studies, comparisons of 
absolute abundances are difficult between different sources using literature values. Assuming the same source sizes 
of CH3OH and complex molecules toward protostars observed with single-dish telescopes result in abundances of a 
few % with respect to CH3OH, comparable to those observed toward Bl-b (Maret ct al. 2005; Bottinelli et al. 2004, 
2007). Where beam averaged or resolved values of complex molecule abundances are actually reported, e.g. toward 
the low-mass outflow L1157 and the hot cores (or accretion disks) of the two IRAS 16293-2422 cores, the abundances 
are also of the order of a few % (Bisschop et al. 2008; Arce et al. 2008). There is thus no evidence that the complex 
molecule abundances with respect to CH3OH toward Bl-b are special compared to other low-mass sources, ranging 
from outflows to hot cores. 

Table ?? lists the abundance ratios of different complex molecules toward the different low-mass sources. The 
Bl-b CH3CHO/HCOOCH3 ratio falls in between what has been previously measured toward the A and B cores in 
IRAS 16293. The upper limits on the CH3OCH3/ HCOOCH3 and C 2 H 5 OH/HCOOCH 3 ratios arc comparable to 
the single-dish observations toward NGC1333 2A, 4A and 4B, and possibly toward the IRAS 16293 envelope as well, 
but significantly lower compared to the resolved observations toward the IRAS 16293 A and B cores. In resolved 
observations toward NGC 1333 2A there is also a CH3OCH3 detection pointing to a higher CH3OCH3 to HCOOCH3 
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ratio in the warm region close to the protostar (J0rgensen et al. 2005). This is entirely consistent with a scenario 
where complex molecules form in cold CO:CH30H ice mixtures in the outer envelope, followed by further formation 
of CH3OHCH3 and C2H5OH in the warmer CO-dcpletcd CH3OH ice close to the protostar (Fig. 7). In the cold parts 
complex ices are only released non-thermally through e.g. UV photodcsorption, while closer to the protostar thermal 
desorption is likely to dominate. This agrees with the low excitation temperatures in the single dish observations 
compared to the observations that resolve the innermost part of the envelope. 

5. CONCLUSIONS 

1. Cold CH3OH gas (excitation temperature of ~10 K) is abundant and widespread toward Bl-b. Quiescent 
CH3OH is most abundant in between the Bl-b SCUBA dust core and the Bl-b protostar, indicative of UV 
photodcsorption of ice on the side of the quiescent dust core because of radiation from the protostar escaping 
through outflow cavities. 

2. Cold HCOOCH3 and CH3CHO are both detected and CH3OCH3 is tentatively detected toward the quiescent 
CH3OH peak. No complex molecules are observed toward an outflow associated with Bl-b, but 3-u upper limits 
with respect to CH3OH are only slightly lower compared to the quiescent core. In addition an asymmetry in the 
CH3CHO emission lines is suggestive of some contribution from a small outflow included in the beam. 

3. Assuming a 10 K excitation temperature the calculated beam-averaged abundances are 2.3% for HCOOCH3, 
1.1% for CH3CHO, <0.8% for CH3OCH3 and <1.0-1.1% for C 2 H 5 OH and HOCH 2 CHO with respect to CH 3 OH. 

4. Building on recent experiments, the observations of large abundances of HCO-containing molecules are explained 
by UV/cosmic ray processing of cold CH30H:CO ice followed by non-thermal desorption of a fraction of the 
produced organic ice. 

5. The beam-averaged abundances with respect to CH3OH and the ratios between different complex molecules are 
similar to other single-dish observations toward low-mass protostars. In contrast resolved observations of the 
cores around protostars are more abundandant in complex molecules without a HCO-group. This is consistent 
with a complex ice chemistry that evolves from HCO-rich to HCO-poor as the ice warms up and the CO-ice 
evaporates when the icy dust falls toward the protostar. 

Including Bl-b, there are still only a handful of observations of complex molecules toward low-mass protostars. 
In addition few complex molecules are detected toward each source. This makes it difficult both to determine the 
prevalence of a complex organic chemistry around low-mass protostars and to firmly establish the main formation 
path of such molecules. The detections of complex molecules in the vicinity of the CHaOH-ice-rich protostar suggests 
that targeting other CHsOH-ice rich protostars may increase our sample of complex molecule detections significantly. 
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Maybe more important from a formation pathway point of view is to increase our understanding on how the complex 
chemistry varies between warm, luke-warm and cold regions i.e. how it varies with distance from the protostar. 
This should be pursued both on larger scales with single dish observations and on smaller scales with interferometers - 
similarly to what has been done toward IRAS 16293 and NGC1333 2a (e.g. J0rgensen et al. 2005; Bisschop et al. 2008). 
In the meantime the complex molecules observed toward Bl-b, and their exceptionally low excitation temperatures 
and line widths, provides clear evidence, for the efficient formation of complex ices during star- and planet-formation. 

We are grateful to the IRAM staff for help with the observations and reduction of the resulting data. This work 
has benefitted from discussions with Herma Cuppen, Robin Garrod and Lars Kristensen and from comments by an 
anonymous referee. Umut Yildiz carried out the second round of observations at the IRAM 30-m with the financial 
support of RadioNct. Support for K. I. is provided by NASA through Hubble Fellowship grant awarded by the 
Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., for NASA, under contract NAS 5-26555. Astrochemistry in Leiden is supported by a SPINOZA grant of the 
Netherlands Organization for Scientific Research (NWO). The Centre for Star and Planet Formation is funded by the 
Danish National Research Foundation and the University of Copenhagen's programme of excellence. 
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Fig. 1. — An integrated CH30H-A + 5os— 4oi map toward the Bl-b region (contours) plotted together with a SCUBA 850 (im map (gray 
scale). The 'P' marks the 'protostar' position, the 'C the 'core' position and 'O' the 'outflow' position. The CH3OH contours are ~2<r i.e. 
0.6 K km s _1 . The map also shows the position of nearby protostar Bl-d and the main outflow knot between the two protostars (diamond) 
as well as possible outflow trajectories from either Bl-b or Bl-d to the knot. 
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Fig. 2. — A channel map of the CH30H-A + 5o5— 4oi line, centered on the Bl-b SCUBA dust core. The contours are ~lcr i.e. 0.1 K and 
the Vi sr of each channel is listed in the upper left corners. The complex molecule observations were pointed toward the CH3OH core (C) 
seen clearest around km s _1 and the south-west outflow (O) visible at <-0.55 km s _1 . The second perpendicular outflow is seen in the 
>1.0 km s _1 channels and is labeled with 02. 
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Fig. 6. — Reaction scheme of a CHsOH-based ice photochemistry adapted from Oberg et al. (2009b) to high-light the impact of CO. In 
a CO-dominated ice the production of several HCO-bearing complex species (red dashed boxes) is enhanced relative to molecules such as 
C2H5OH, CH3OCH3 and (CH20H)2, which tend to dominate among the pure CHsOH-photochcmistry products (solid and dotted boxes). 
(CH20H)2 and HOCH2CHO (dotted boxes) are only produced abundantly at luke-warm temperatures (>30 K in the laboratory). 
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TABLE 1 
Observational parameters 



Positions and CH3OH observations with HERA 




Name a 


R.A. (J2000) Dec. (J2000) rms (mK) 8V 


(km s 


- 1 ) 




'core' 


03:33:20.80 31:07:40.0 


92 


0.39 






'outflow' 


03:33:21.90 31:07:22.0 


65 


0.77 






'protostar' 


03:33:20.34 31:07:21.4 


92 


0.39 




Complex molecules observations with EMIR 


Setting 


Receiver 


Frequency band (GHz) rms 


(mK) 


8V 


(km s 1 ) 




'core' 


'outflow' 


'core' 


'outflow' 


1 


E90 


88.81-88.93 5.4 


2.0 


0.26 


2.1 




E150 


146.84-146.96 5.5 


2.2 


0.32 


1.3 


2 


E90 


115.46-115.58 14.6 


9.9 


0.41 


1.6 




E150 


146.93-147.05 4.6 


3.9 


0.32 


1.3 


3 


E150 


130.82-130.94 8.4 




0.18 





a The 'core' is the quiescent CH3OH maximum ~10 offset from the SCUBA core and the 'protostar' is the source discovered by Spitzer. 



TABLE 2 
Observed CH3OH lines with 1-ct uncertainties in brackets. 



Position 


Transition 


Freqency 


Slow 


FWHM 


fTdV 10" 


fTdV 19" 


fTdV 27" 








GHz 


cm 


km s _1 


K km s 


K km s 


K km s — 1 


'core' 


A+5 05 


-4 4 


241.7914 


16.1 


1.2[0.1] 


2.8[0.6] 


1.5[0.3] 


1.2[0.3] 




E-5i 8 - 


-4i4 


241.7672 


20.0 


1.3[0.1] 


2.2[0.4] 


1.2[0.2] 


1.1[0.2] 




E-5o5 - 


-4 04 


241.7002 


25.3 


2.4[0.5] 


0.5[0.1] 


0.2[0.1] 


0.2[0.1] 


'outflow' 


A+5,)5 


-4 4 


241.7914 


16.1 


4.0[0.2] 


3.3[0.7] 


1.9[0.4] 


1.4[0.4] 




E-5i5 - 


-4i4 


241.7672 


20.0 


4.4[0.3] 


2.8[0.5] 


1.6[0.3] 


1.2[0.3] 




E-5 05 - 


-4 4 


241.7002 


25.3 


5.7[1.6] 


0.5[0.1] 


0.2[0.1] 


0.2[0.1] 


'protostar' 


A+5o5 - 


-4 4 


241.7914 


16.1 


1.1[0.1] 


1.6[0.3] 


1.2[0.3] 


1.2[0.3] 




E-5i5- 


-4i4 


241.7672 


20.0 


1.0[0.1] 


1.4[0.3] 


1.1[0.2] 


1.0[0.2] 




E"5o5 - 


-4 04 


241.7002 


25.3 


1.0[0.9] 


0.5[0.1] 


0.2[0.1] 


0.2[0.1] 
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TABLE 3 
Observed complex molecular lines with 1-ct uncertainties in brackets. 



Molecule 


Transition 


Frcq. 


Blow 


FWHM a 


fTdV b 








GHz 


cm -1 


km s -1 


mK km s -1 


HCOOCH3-E 


7i 6 - 


- 7i 5 


88.8432 


9.5 


0.55[0.13] 


18 [4] / < 24 


HCOOCH3-A 


7l6" 


- 7i 5 


88.8948 


9.5 


0.60[0.12] 


16 [4] / < 24 


C2H5CN 


19l 18 - 


- 19(, 19 


88.8948 


55.6 




< 11 / < 24 


HOCH2CHO 


9 4 6~ 


-9 s7 


88.8924 


21.2 




< 11 / < 24 


CH3CHO-A 


62 5 - 


-52 4 


115.4939 


15.9 




< 30 / < 119 


C2H5OH 


515- 


-50 5 


115.4945 


48.2 




< 30 / < 119 


CH3OCH3 


5l 5 2 - 


- 5(J4 2 


115.5440 


6.3 




< 30 / < 119 


CH3OCH3 


5l 5 3 - 


- 50 4 3 


115.5440 


6.3 




< 30 / < 119 


CH3OCH3 


5i 5 1 - 


- 50 4 1 


115.5448 


6.3 




< 30 / < 119 


CH3OCH3 


5i 5 - 


- 5o 4 


115.5457 


6.3 




< 30 / < 119 


C2H5CN 


633 - 


-624 


115.5474 


9.4 




< 30 / < 119 


C2H5OH 


43 2 2 - 


-4 2 32 


130.8715 


9.3 




< 17 


CH3CHO-E 


7l7" 


- 61 6 


130.8918 


14.7 


0.89[0.18] 


39 [6] 


CH3CHO-A 


7it- 


- 61 6 


130.8927 


14.7 


0.67[0.18] 


42 [6] 


C2H5CN 


15o 15 - 


~ 14() 14 


130.9039 


30.9 




< 17 


CH3OCH3 


53 32 - 


"5242 


146.8660 


13.4 




< 11 / < 26 


CH3OCH3 


53 33 " 


-5243 


146.8703 


13.4 




< 11 / < 26 


CH3OCH3 


5331 " 


- 52 4 1 


146.8725 


13.4 




< 11 / <26 


CH3OCH3 


53 30 " 


"5240 


146.8773 


13.4 




< 11 / < 26 


NH2CHO 


7o 7 8 " 


- 60 6 7 


146.8716 


14.8 




< 11 / < 26 


NH2CHO 


7o 76 " 


- 60 6 5 


146.8716 


14.8 




< 11 / < 26 


NH2CHO 


7o 7 7 " 


- 60 6 6 


146.8717 


14.8 




< 11 / < 26 


C2H5CN 


17i 17 - 


- 16l 16 


146.8826 


40.1 




< 11 / < 26 


HCOOCH3-E 


12s 10 


-113 9 


146.9777 


31.3 




< 9 / < 47 


HCOOCH3-A 


123 10 


-H39 


146.9880 


31.2 




< 9 / < 47 


CH3OCH3 


7l 7 2 ' 


- 60 6 2 


147.0242 


13.2 




< 9 / < 47 


CH3OCH3 


7i 7 3 - 


- 60 6 3 


147.0242 


13.2 




< 9 / < 47 


CH3OCH3 


7l 7 1 - 


- 60 6 1 


147.0249 


13.2 


~0.53 


~8 / < 47 


CH3OCH3 


7l 70 - 


- 60 6 


147.0256 


13.2 




< 9 / < 47 



1) 



FWHM for detected lines toward the 'core' position. 
b Values arc for 'core' / 'outflow' positions. 3-ct upper limits are calculated from the rms assuming line widths of 0.68 (average of the four 

complex line detections) and 4.2 km s for the 'core' and 'outflow' positions, respectively. 
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TABLE 4 
Calculated abundances and 3-ct upper limits with l-o- uncertainties in brackets. 



Position 


Species 


T rot / K 


N a / 10 13 cnr 2 


% CH 3 OH a 




core' 


CH3OH 


10[5] 


~36 / 47 b 


100 




core' 


HCOOCH3-E 


- 


0.44[0.10] 


1.2[0.3] 




core' 


HCOOCH3-A 


- 


0.39(0.10] 


1.1[0.3] 




core' 


CH3CHO-E 


- 


0.26[0.04] 


0.6[0.1] 




core' 


CH3CHO-A 


- 


0.28[0.04] 


0.6[0.1] 




core' 


CH3OCH3 


- 


5-0.36 


$0.8 




core' 


C2H5OH 


- 


<0.48 


<1.0 




core' 


HCOCH2OH 


- 


<0.41 


<1.1 




outflow' 


CH3OH 


8[4] 


~78 / 83 b 


100 




outflow' 


HCOOCH3-E 


- 


<0.56 


<0.7 




outflow' 


HCOOCH3-A 


- 


<0.56 


<0.7 




protostar' 


CH3OH 


11[1] 


23 [4] 


100 



a Assuming an excitation temperature of 10 K for the complex species toward both the core and the outflow position. 
b The first value is derived assuming a 27" beam and the second assuming a 19" beam. 



TABLE 5 
Calculated abundances and 3-<t upper limits with 1-cr uncertainties in brackets. 



Source 



T rot FWHM CH3CHO CH3OCH3 C2H5OH 

K km s- 1 / HCOOCH 3 a / HCOOCH 3 a / HCOOCH 3 a 

<0.4 



2/1 



-b core 


10[10] 


0.55-0.89 


0.5[0.3] 


$0.3 


AS 16293 cnv. b 


14-67 


0.8-10 


0.1 


0.67 


AS 16293 core A / B c 


100-300 


1.5-2.5 


<0.1 / 3 


1/2 


5C1333 IRAS 2A d >° 


38-101 


1 1 


- 


<0.3 


5C1333 IRAS 4A d 


24-36 


1.2-4.2 


- 


<0.4 


^C1333 IRAS 4B d 


34-38 


0.5-1.8 


- 


<0.7 


157 outflow' 


18-27 


2.3-5.5 


_ 


_ 



0.4 



a Whcre only A or E column densities are available the total HCOOCH3 abundance is dervied by assuming equal abundances of both. 
b Cazaux ct al. (2003); Herbst & van Dishoeck (2009) 

c Beamsizc of 5.5" x3.2" Bisschop ct al. (2008), 1.3" x2.7" Kuan ct al. (2004) and Huang ct al. (2005) 
Abundances are rescalcd assuming the same beam dilution for CH3OH and complex organic molecules. Single dish column densities 

from Marct ct al. (2005); Bottinelli ct al. (2004, 2007) 
C CH30CH3 is detected toward the central core using the Sub-Millimeter Array (j0rgensen et al. 2005). 
f Arcc ct al. (2008) 
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Fig. 8. — Observed lines in the five different settings explored toward the 'core' position plotted versus AV, the deviation from the core 
velocity 6.3 km s — 1 . The spectra are centered on the rest frequency shown in each panel. 
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Fig. 9. — Observed lines in the five different settings explored toward the 'outflow' position plotted versus AV, the deviation from the 
core velocity 6.3 km s . The spectra are centered on the rest frequency shown in each panel. 



